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Abstract 

Transition-metal (TM)-doped diluted magnetic oxides (DMOs) have attracted attention from both experimental and 
theoretical points of view due to their potential use in spintronics towards new nanostructured devices and new 
technologies. In the present work, we study the magnetic properties of Sn 0 9 6 TM 004 O2 and Sn 096 Jf^\ 0m O^ 98 {V o ) 0mi 
where TM = Fe and Co, focusing in particular in the role played by the presence of 0 vacancies nearby the TM. The 
calculated total energy as a function of the total magnetic moment per cell shows a magnetic metastability, 
corresponding to a ground state, respectively, with 2 and 1 |u B /cell, for Fe and Co. Two metastable states, with 0 
and 4 |j B /cell were found for Fe, and a single value, 3 Me/cell, for Co. The spin-crossover energies (E s ) were 
calculated. The values are E% /2 = 107 meV and /F 4 / 2 = 25 meV for Fe. For Co, /is 71 = 36 meV. By creating 0 vacancies 
close to the TM site, we show that the metastablity and E s change. For iron, a new state appears, and the state 
with zero magnetic moment disappears. The ground state is 4 |j B /cell instead of 2 |i B /cell, and the energy Ef 4 is 30 
meV. For cobalt, the ground state is then found with 3 |i B /cell and the metastable state with 1 |j B /cell. The 
spin-crossover energy E\ /3 is 21 meV. Our results suggest that these materials may be used in devices for spintronic 
applications that require different magnetization states. 
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Background 

Nowadays, dilute magnetic oxides (DMOs) are potential 
candidates for both spintronic devices and nanodevices 
applications. Although the existence of room temperature 
ferromagnetism (FM) in transition metal (TM) -doped 
Sn0 2 has been reported, the origin of the FM is still con- 
troversial There are indications that the FM comes from 
different sources, metallic clusters, secondary phases, or 
is due to a free carrier-mediated mechanism in the bulk. 
The presence of oxygen vacancies is systematically related 
to the observed ferromagnetic state. These systems are 
good candidates to obtain materials with a half-metallic 
behavior, with 100% spin-polarized carriers at the Fermi 
level. 

Tin dioxide (Sn0 2 ) doped with transition metals 
(TMs) have been extensively investigated recently due to 
the resulting important magnetic properties [1-6]. The 
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ferromagnetic behavior has been observed at room 
temperature in Cr-, Mn-, Fe- and Co-doped Sn0 2 DMO 
systems [7-12], indicating the potential of such systems 
for spintronic applications. It has also been observed 
that the presence of oxygen vacancies appears to be 
required for producing FM in DMOs, such as, e.g., in 
Co-doped ZnO [13], in Co-doped Ti0 2 [14,15], in Fe- 
and Co-doping in ln 2 0 3 [16,17], and in Fe-, Co-, and 
Cr- doped Sn0 2 [7,18,19]. A theoretical model proposed 
by Coey et al. to interpret the FM in these semiconduct- 
ing oxides requires the existence of oxygen vacancies in 
close proximity to TM sites in order to maintain the 
charge neutrality [20]. For Cr-doped Sn0 2 nanoparticles, 
the obtained FM behavior is limited by a maximum dop- 
ing concentration x L which has a strong relation with 
structural changes revealed from X-ray diffraction mea- 
surements [7]. The presence of oxygen vacancies in 
these Sn! _ ^TM^0 2 samples, in which the TM concen- 
trations x varies from 0% to 10%, has been detected by 
electron paramagnetic resonance experiments [19]. 
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Many efforts have been made in attempt to characterize 
and understand the mechanisms involved in the 
ferromagnetic behavior observed in such systems. 

In this work, we study the oxygen vacancy influence 
on magnetic properties of the Fe- and Co-doped Sn0 2 
diluted alloys. First, the systems Sn 2 _ x MT x 0 2 , for x = 
0.04, were studied through ab initio electronic structure 
calculations performed within the spin density functional 
theory. The concentration of x = 0.04 corresponds to a 
typical experimental value. Second, an oxygen vacancy 
nearest neighbor to the TM atom in the alloys was intro- 
duced. Sn x _ x TM x 0 2 - y(Vo) y systems, with x = 0.04 
and y = 0.02, and its consequences for the magnetic 
behavior of these systems were considered. Finally, an 
investigation about the magnetic metastability was 
done, and the spin-crossover phenomenon was studied. 
A metamagnetic state is the key underlying conceptual 
mechanism for storage, memory, and display device 
and nanodevice applications [21], and this work shows 
that DMO materials based on Fe- and Co-doped Sn0 2 
taking into account the oxygen vacancy influence could 
be engineered to display different stable magnetized 
states. 

Methods 

The calculations were based on the spin density 
functional theory. We employed the projector augmen- 
ted wave method implemented in the Vienna ab-initio 
simulation package (VASP-PAW) [22,23]. The exchange- 
correlation potential used was the generalized gradient 
approximation in the Perdew, Burke, and Ernzerhof 
[24] approach. The method has been previously used to 
study the structural and electronic properties of bulk 
rutile Sn0 2 [25]. The valence electronic distributions 
for the PAWs representing the atoms were Sn 
4d 10 5s 2 5/, Fe 3rf 7 4s\ Co 3d 7 4s 2 , and O 2s 2 2p 6 . The 
onsite correction Hubbard U for Co- and Fe d orbitals 
was not considered. Previous calculations for Cr doping 
Sn0 2 taking into account a U correction show no 
major changes in our conclusions. Scalar relativistic 
effects were taken into account. To describe the alloys, 
we used a 72-atom supercell (24 Sn and 48 O atoms) 
and a 4 x 4 x 4 mesh of Monkhorst-Pack k-points for 
integration in the Brillouin zone. All the calculations 
were done with a 490-eV energy cutoff in the plane- 
wave expansions, and the systems were fully relaxed 
until the residual forces on the ions were less than 10 
meV/A. 

Results and discussion 

We studied the systems Sn 0 96 TM 0 . 04 O 2 and the 
Sn 0 .96TMo.o40i. 98 (Vo)o.o2 (with TM = Fe and Co) with 
the oxygen vacancy as the TM nearest neighbor. In both 
cases, a single tin atom was substituted with a TM atom 



in the 72-atom supercell, simulating the x = 0.04 impur- 
ity and y = 0.02 oxygen vacancy concentrations. For all 
cases, the total energy was calculated for several mag- 
netic moment values. 

It was recently shown by us [26] that when a Sn atom 
is replaced by a chromium atom in Sn0 2 , a high-spin 
(HS) ground state with a magnetic moment m = 2 u B / 
cell and a low-spin state (LS) with a magnetic moment 
m = 0 [i B /cell are obtained. For this case, a spin cross- 
over becomes possible with an energy barrier of 114 
meV calculated for the transition from m = 0 to 2 \i B / 
cell. When an oxygen vacancy (out of the six first 
neighbors of Cr) is considered, the behavior of the total 
energy versus the magnetic moment per cell showed 
the appearance of a second HS configuration, with 
magnetic moment m = 4 [i B /cell and an energy barrier 
of 32 meV relative to the 2 [i B /cell state. The ground 
state, however, remains as the 2 |i B /cell magnetic mo- 
ment HS state. The energy barrier for the m = 0 to 2 
(i B /cell transition was reduced to 27 meV. Comparing 
this value with the one for this barrier in the DMO 
without the vacancy (114 meV), a drastic transition for 
the DMO without the vacancy of 114 meV, a drastic re- 
duction by about 75% is observed. 

We show here that if a single Sn atom is substituted 
by a Fe or a Co impurity, the magnetic metastability 
observed for Cr-doped SnO$_2$ diluted alloys also 
occurs. Figure la,b shows the total energy behavior for 
the Sn 0 .96Fe 0 .04O 2 and Sno.96Fe 0 .o40i.98(Vb)o.o2 diluted 
magnetic alloys, respectively, in function of magnetic 
moment/cell. A ground state with a magnetic moment 
m = 2 (i B /cell and two metastable states with magnetic 
moments of m = 0 and m = 4 |i B /cell were observed for 
the system without the oxygen vacancy. The crossover 
energy barrier for the transition from m = 0 to 2 |i B /cell 
is £^ /2 = 107 meV, and for the transition from m = 4 to 
2 ne/cell is £f/ 2 = 25 meV. When the oxygen vacancy is 
considered, m = 0 |i B /cell is no longer a metastable state; 
however, new states appear with magnetic moment vary- 
ing in the range 2 < m < 6. The ground state occurs at 
m = 4 (i B /cell, and two almost flat regions between 2 < 
m < 3.5 (i B /cell (region Rl) and 4 < m < 6 [i B /cell (region 
R2) are seen. Energy barriers of £?/ 4 = 30 meV for the 
transition from m = 2 to 4 (i B /cell and £| /4 =16 meV for 
the transition from m = 4 to 6 |i B /cell were observed. 
Figure lc,d shows the total energy versus magnetic mo- 
ment for Sn 0 .96Coo.o40 2 and Sn 0 .96Coo.o40 1 .98(V 0 )o.o2 
alloys, respectively. The ground state for the system 
without the oxygen vacancy is at m = 1 |i B /cell, and a 
magnetic metastable at m = 3 |i B /cell is seen. The energy 
barrier of E^ 1 = 36 meV for the transition from m - 3 to 
1 (i B /cell was obtained. For the Co-doped Sn0 2 with an 
oxygen vacancy, we obtain a change in the magnetic 
metastable states: the ground state is at m = 3 |i B /cell 
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Figure 1 Total energy vs. magnetic moment per cell. For (a) Sn 0 .96Feo.o40 2 , (b) Sn 0 .96Fe 0 .o40 1 .9 8 (i/o)o.o2/ ( c ) Sno.9eCoo.04O2, and (d) 
Sn 0 .96Co 0 .o40 1 .9 8 (Vb)o.o2. The total energy of the high-spin ground state is set to zero. 



and a magnetic metastable at m = 1 |i B /cell is seen, with 
a crossover barrier energy of El, 13 = 21 meV from m = 1 
to 3 |i B /cell. 

A better understanding of the behavior for the mag- 
netic moment can be obtained if we analyze the oxida- 
tion states of the TM and of the neighbor atoms around 
it. For the Sn 0 .96Fe 0 .04O 2 alloy, the Fe 4+ (3d 4 ) impurity 
replacing Sn 4+ allows the magnetic states m = 0, 2, and 
4 (i B /cell. In this case, the neighborhood does not con- 
tribute to the magnetization of the system. When an 
oxygen vacancy is taken into account, the impurity state 
changes from Fe 4+ to Fe 3+ (3d 5 ) and to Fe 2+ (3d 6 ), and 
the neighboring atoms contribute to the magnetic mo- 
ment. For m = 2 |i B /cell, the oxidation state is Fe 3+ 
where four spin-up and one spin-down electrons from 
Fe plus one spin-down electron arising from the neigh- 
boring atoms allow for the metastable state (with m = 2 
(i B /cell). For the m = 4 and 6 |i B /cell, the oxidation state 
is Fe 2+ where five spin-up and one spin-down electrons 
from Fe originate the state m = 4 |i B /cell, while five 
spin-up and one spin-down electrons from Fe plus two 
spin-up electrons arising from the neighboring atoms 
allow the metastable state with m = 6 |i B /cell. These 
findings are in agreement with the experimental data 
[27-30]. 



Likewise, for the Sn 0 96 Co 0 . 04 O2 dilute magnetic alloy, 
the Co 4+ (3d 5 ) and Co 3+ (3d 6 ) impurity replacing Sn 4+ 
give rise to the magnetic states m = 1 and 3 |i B /cell, re- 
spectively. For m = 1 |i B /cell, three spin-up and two 
spin-down electrons from cobalt and for m = 3 |i B /cell 
four spin-up and two spin-down electrons from cobalt 
plus one spin-up electron arising from the neighboring 
atoms allow this metastable state. The states Co 3+ (3d 6 ) 
and Co 2+ (3d 7 ) are possible when an oxygen vacancy is 
taken into account for m = 1 and 3 |i B /cell, respectively. 
For m = 1 |i B /cell, four spin-up and two spin-down elec- 
trons plus one spin-down electron from neighboring 
atoms are involved. For m = 3 (_i B /cell, five spin-up and 
two spin-down electrons from cobalt give rise to this 
metastable state. Experimental studies have confirmed 
the incorporation of Co 2+ cations into the rutile Sn0 2 
lattice [31]. 

As observed previously for the chromium impurity 
[26], for iron and cobalt we also observe a relationship 
between the structural modification around the TM 
atom, due to the electronic and ionic relaxations, and 
the occurrence of the magnetic metastability. If we con- 
sider a spherical volume involving the TM whose radius 
is an average distance between the TM and the six oxy- 
gen nearest neighbors, our calculations showed that, 
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after full relaxation, the corresponding volume is 
reduced. For Co, the volume is reduced by about 16% 
for the HS state (2 u B /cell) and by 17.5% for the LS state 
(0 u B /cell). Considering the oxygen vacancy, for the HS 
states (2 and 4 u B /cell) the volume reductions, after full 
relaxed calculations, were 20% and 22%, respectively, 
while for the LS state (0 u B /cell), it was 31%. For LS con- 
figurations, the presence of an oxygen vacancy allows 
greater relaxations which reduce the total energy of the 
system lowering the energy barrier for the crossover. As 
shown in Figure 2a for Fe, the volume is reduced by 
about 19% for the states m = 0 and 2 u B /cell, and by 
14% for the state m = 4 u B /cell. If the oxygen vacancy is 
present, the volume around the Fe impurity is reduced 
by about 19% for the magnetic moment between 2 and 
3.5 u B /cell, region Rl, and by 12% for the range between 
4 and 6 u B /cell, region R2, shown in Figure 2b, both cor- 
responding to an almost flat region. For the cobalt 



a 





91 -r 




90- 




89- 




88- 




87- 








86- 


CD 




E 


85- 


=5 


84- 


O 


> 


83- 




82- 




81- 



Sn 0. 9 6 Fe 0.04 O 2 



"1 1 1 - 



1 2 3 4 5 

Magnetic Moment (n /cell) 



89 
88 
87 
86 
85 
84 
83 

3 81 
80 



CD 



Sn 0,6 Fe 0.04 O 1.98( V o)0. 



R2 



R1 



79 
78 



0 1 2 3 4 5 6 

Magnetic Moment (jn /cell) 

Figure 2 Relative volume (%) around the iron atom vs. 
magnetic moment per cell. For (a) Sn 0 .96Feo.o40 2 and (b) 
Sno.96Feo.04O1. 98(^0)0.02- 



impurity, the volume around Co is reduced by about 
19% for the state m = 1 u B /cell and by 15% for the state 
m = 3 u B /cell. Considering an oxygen vacancy near the 
Co atom, the volume is reduced by about 19% for the 
state m = 1 u B /cell and by 11% for the state m = 3 u B / 
cell. The observed magnetic metastability in the DMOs 
studied here is attributed to a structural modification 
(relaxations) around the TM impurity. The presence or 
absence of magnetism in this case is determined by a 
competition between intra-atomic exchange interactions 
and inter- atomic electronic motion due to the crystalline 
field [32]. Therefore, the occurrence of the LS and HS 
states depends on the effective balance between these 
two interaction fields. The LS state occurs when the 
crystal field splitting is larger than the intra-atomic ex- 
change splitting (lowest volume); otherwise, the HS state 
is the ground state (highest volume). 

The ground state total and projected density of states of 
the Fe 3d orbital are shown in Figure 3a,b for 
Sno.9eFeo.04O2 and Sn 0 .96Feo.o40i.9 8 (Vo)o.o2 alloys, respect- 
ively. The vertical lines represent the energy values of the 
highest occupied states for the majority and minority 
spins. For both cases, the iron 3d-derived states are found 
to lie in the gap region, showing a half-metallic behavior. 
For Sno.96Coo.04O2 and Sn 0 .96Coo.o40i. 9 8(Vo)o.o2 alloys, the 
DOS are shown in Figure 3c,d, respectively. For all cases, 
we observe a strong character arising from the 3d orbital 
near the gap region. This feature is also responsible for 
the half-metallic behavior displayed by the DMOs studied 
here. 

Conclusions 

The influence of the oxygen vacancy in the magnetic 
and electronic properties of iron and cobalt as impur- 
ities in a DMO configuration in rutile Sn0 2 was studied 
using first principle calculations performed within the 
spin-density functional theory. A magnetic metastability 
was observed for both impurity cases. Energy barriers 
were obtained for the spin crossover between the m = 
0 and 2 u B /cell and between the m = 2 and 4 u B /cell 
states in Sn 0 .96Fe 0 .04O 2 . For Sn 0 .96Co 0 .04O 2 , the observed 
magnetic metastability and energy barrier were obtained 
for the spin crossover between the m = 3 and 1 u B /cell 
states. 

When an oxygen vacancy is considered as one of the six 
first neighbors to the Fe and Co impurities in these alloys 
(Sn 0 .96Feo.o40i.9 8 (Vo)o.o2 and Sn 0 .96Coo.o40i.9 8 (Vo)o.o2) ) a 
considerable modification is observed in the magnetic 
metastability behavior, with new allowed states appear- 
ing for iron. For the cobalt impurity, the presence of 
an oxygen vacancy promoted a ground state exchange 
between the magnetic state 1 and 3 u B /cell. This 
behavior is attributed to the relative contributions of 
the intra-atomic exchange interaction effects and the 
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Figure 3 The ground state total (black lines) and projected (gray shaded areas) density of states. The transition metal 3d-derived states 
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inter-atomic electron motion effects due to the crystal- 
line field, which are responsible for the relaxations 
around the TM impurities. Finally, to manipulate the 
electron spin opens new possibilities to engineering new 
spintronic devices, and TM-doped DMOs are potential 
candidates to represent a new kind of material, which 
display magnetic metastability, SCO phenomena, and a 
half-metallic behavior. 
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